AFRL-RB-WP-TP-2  011-3104 


ALTITUDE  CONTROL  OF  A  SINGLE  DEGREE  OF 
FREEDOM  FLAPPING  WING  MICRO  AIR  VEHICLE 
(POSTPRINT) 

David  B,  Doman,  Michael  W,  Oppenheimer,  Michael  A,  Bolender,  and  David  O,  Sigthorsson 

Control  Design  and  Analysis  Branch 
Control  Sciences  Division 


AUGUST  2009 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 


STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
AIR  VEHICLES  DIRECTORATE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7542 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a 
currently  valid  0MB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YY) 

August  2009 

2.  REPORT  TYPE 

Conference  Paper  Postprint 

3.  DATES  COVERED  (From  ■  To) 

05  November  2008  -  13  August  2009 

4.  TITLE  AND  SUBTITLE 

ALTITUDE  CONTROL  OF  A  SINGLE  DEGREE  OF  FREEDOM  FLAPPING  WING 
MICRO  AIR  VEHICLE  (POSTPRINT) 

5a.  CONTRACT  NUMBER 

In-house 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6220 IF 

6.  AUTHOR(S) 

David  B.  Doman,  Michael  W.  Oppenheimer,  Michael  A.  Bolender, 
and  David  0.  Sigthorsson  (AFRL/RBCA) 

5d.  PROJECT  NUMBER 

2401 

5e.  TASK  NUMBER 

N/A 

5f.  WORK  UNIT  NUMBER 

QI2K 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Control  Design  and  Analysis  Branch  (AFRL/RBCA) 

Control  Sciences  Division 

Air  Force  Research  Laboratory,  Air  Vehicles  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7542 

Air  Force  Materiel  Command,  United  States  Air  Force 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-RB-WP-TP-20 1 1  -3 1 04 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 

Air  Vehicles  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7742 

Air  Force  Materiel  Command 

United  States  Air  Force 

10.  SPONSORING/MONITORING 
AGENCY  ACRONYM(S) 

AFRL/RBSD 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-RB-WP-TP-20 11-3104 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

13.  SUPPLEMENTARY  NOTES 

PAO  Case  Number:  88ABW-2009-3332;  Clearance  Date:  21  Jul  2009.  Document  contains  color. 

Conference  paper  published  in  the  proceedings  of  the  AIAA  Guidance,  Navigation,  and  Control  Conference  held  10-13 
August  2009  in  Chicago,  IL. 

14.  ABSTRACT 

A  control  strategy  is  proposed  for  a  minimally  actuated  flapping  wing  micro  air  vehicle.  The  Harvard  RoboFly  vehicle 
accomplished  the  first  takeoff  of  an  insect  scale  flapping  wing  aircraft.  This  flight  demonstrated  the  capability  of  the 
aircraft  to  accelerate  vertically  while  being  constrained  by  guide-wires  to  avoid  translation  and  rotation  in  the  other  five 
degrees  of  freedom.  The  present  work  proposes  an  altitude  control  scheme  that  would  enable  a  similar  vehicle  under  the 
same  constraints  to  hover  and  track  altitude  commands.  Using  a  blade  element-based  aerodynamic  model  and  cycle 
averaging,  it  will  be  shown  that  altitude  control  of  such  an  aircraft  can  be  achieved.  The  RoboFly  makes  use  of  a  single 
bimorph  piezoelectric  actuator  that  symmetrically  varies  the  angular  displacement  of  the  left  and  right  wings  in  the  stroke 
plane.  The  wing  angle-of-attack  variation  is  passive  and  is  a  function  of  the  instantaneous  angular  velocity  of  the  wing  in 
the  stroke  plane.  The  control  law  is  designed  to  vary  the  frequency  of  the  wing  beat  oscillations  to  control  the  longitudinal 
body-axis  force  which  is  used  to  achieve  force  equilibrium  in  hover  and  acceleration  when  tracking  time-varying  altitude 
commands. 


15.  SUBJECT  TERMS 

flapping  wing  micro  air  vehicles,  MAV,  minimal  actuation,  altitude  control 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT: 

SAR 

18.  NUMBER 

OF  PAGES 

26 

19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

ILt  Zachary  H.  Goff 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

N/A 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


AIAA  Guidance,  Navigation,  and  Controi  Conference 
10-13  August  2009,  Chicago,  iilinois 


AIAA  2009-6159 


Altitude  Control  of  a  Single  Degree  of  Freedom  Flapping  Wing 

Micro  Air  Vehicle 


David  B.  Doman  *  Michael  W.  Oppenheimer  |  Michael  A.  Bolender  land  David  O.  Sigthorsson  ^ 


A  control  strategy  is  proposed  for  a  minimally  actuated  flapping  wing 
micro  air  vehicle.  The  Harvard  RoboFly  vehicle  accomplished  the  first 
takeoff  of  an  insect  scale  flapping  wing  aircraft.  This  flight  demonstrated 
the  capability  of  the  aircraft  to  accelerate  vertically  while  being  constrained 
by  guide-wires  to  avoid  translation  and  rotation  in  the  other  five  degrees 
of  freedom.  The  present  work  proposes  an  altitude  control  scheme  that 
would  enable  a  similar  vehicle  under  the  same  constraints  to  hover  and 
track  altitude  commands.  Using  a  blade  element-based  aerodynamic  model 
and  cycle  averaging,  it  will  be  shown  that  altitude  control  of  such  an  aircraft 
can  be  achieved.  The  RoboFly  makes  use  of  a  single  bimorph  piezoelectric 
actuator  that  symmetrically  varies  the  angular  displacement  of  the  left  and 
right  wings  in  the  stroke  plane.  The  wing  angle-of-attack  variation  is  passive 
and  is  a  function  of  the  instantaneous  angular  velocity  of  the  wing  in  the 
stroke  plane.  The  control  law  is  designed  to  vary  the  frequency  of  the 
wing  beat  oscillations  to  control  the  longitudinal  body-axis  force  which  is 
used  to  achieve  force  equilibrium  in  hover  and  acceleration  when  tracking 
time-varying  altitude  commands. 


*Senior  Aerospace  Engineer,  Control  Design  and  Analysis  Branch,  2210  Eighth  Street,  Ste.  21,  Air  Force 
Research  Laboratory,  WPAFB,  OH  45433-7531  Email  David.Doman@wpafb.af.mil,  Ph.  (937)  255-8451,  Fax 
(937)  656-4000,  Associate  Fellow  AIAA 

^Senior  Electronics  Engineer,  Control  Design  and  Analysis  Branch,  2210  Eighth  Street,  Ste  21,  Air  Force 
Research  Laboratory,  WPAFB,  OH  45433-7531  Email  Michael.Oppenheimer@wpafb.af.mil,  Ph.  (937)  255- 
8490,  Fax  (937)  656-4000,  Senior  Member  AIAA 

1  Aerospace  Engineer,  Control  Design  and  Analysis  Branch,  2210  Eighth  Street,  Ste  21,  Air  Force  Research 
Laboratory,  WPAFB,  OH  45433-7531  Email  michael.bolender@wpafb.af.mil,  Ph.  (937)  255-8492,  Fax  (937) 
656-4000,  Senior  Member  AIAA 

^Electronics  Engineer,  Control  Design  and  Analysis  Branch,  2210  Eighth  Street,  Ste  21,  Air  Force  Re¬ 
search  Laboratory,  WPAFB,  OH  45433-7531  Email  David.Sigthorsson@afmcx.net,  Ph.  (937)  255-9707,  Fax 
(937)  656-4000.  This  research  was  performed  while  this  author  held  a  National  Research  Council  Research 
Associateship  Award  at  the  Air  Force  Research  Laboratory. 


Approved  for  public  release;  distribution  unlimited. 

1  of  19 


This  material  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject  to  copyright  protection  in  the  United  States. 


I.  Introduction 


The  first  takeoff  of  an  insect  scale  biomimetic  flapping  wing  micro  air  vehicle  was  achieved 
by  an  aircraft  called  RoboFly  that  was  developed  at  Harvard  University  by  Wood  et.ald  A 
key  feature  that  led  to  the  successful  hrst  flight  is  that  the  vehicle  is  minimally  actuated 
and  makes  use  of  passive  wing  rotation  to  mimic  the  wing  beat  patterns  of  a  dipterian 
insect.  As  shown  in  Figure  1,  RoboFly  uses  a  single  bimorph  piezoelectric  actuator  to 
impart  symmetric  motion  to  two  wings  simultaneously.  Tangential  motion  of  the  tip  of  the 
piezoelectric  actuator  is  converted  to  rotational  motion  of  the  wings  by  way  of  a  linkage.  The 
linkage  elements  are  designed  to  achieve  impedance  matching  between  the  wing  and  actuator 
forces  and  to  amplify  the  relatively  small  motion  of  the  tip  of  the  bimorph  strip  into  large 
angular  displacements  of  the  wing  in  the  stroke  plane.  A  wing  is  connected  to  the  movable 
wing  root  by  a  flexible  hinge  that  provides  for  passive  rotation  of  the  wing.  The  hinge  allows 
the  wing  to  passively  invert  its  orientation  about  the  hinge  joint  as  the  wing  reverses  direction 
at  the  end  of  each  stroke.  This  results  in  planform  rotations  that  approximate  the  motion 
of  dipterian  insect  wings.  As  the  wing  traverses  through  the  stroke  plane,  dynamic  pressure 
acting  on  the  wing  tends  to  cause  it  feather  into  the  wind;  however,  as  shown  in  Figure  2, 
the  hinge  joint  is  designed  for  interference  between  the  planform  and  root  to  prevent  the 
wing  from  over-rotating.  Therefore,  the  wing  holds  a  constant  angle  of  attack  relative  to  the 
stroke  plane  once  a  critical  dynamic  pressure  is  reached.  The  actuator  and  the  carbon  hber 
substrate  to  which  it  is  mounted  are  cantilevered  to  the  fuselage  and  together  with  the  wing 
form  a  spring-mass-damper  system  that  has  a  known  resonant  frequency.  In  the  Harvard 
experiment,  this  dynamic  system  was  driven  at  resonance  for  maximum  energy  efficiency 
to  achieve  flight.  The  hrst  hight  resulted  in  unregulated  hight  up  a  wire  that  constrained 
the  vehicle  motion  to  vertical  translation.  In  the  hrst  part  of  this  paper,  we  investigate  the 
suitability  of  wingbeat  frequency  modulation  to  allow  tracking  of  a  desired  vertical  position 
prohle  of  the  single  degree  of  freedom  (DOF)  RoboFly  experiment.  In  the  companion  pair 
of  papers,^’ ^  a  vehicle  concept  and  associated  control  strategies  that  will  allow  the  vehicle 
to  break  free  of  the  wire  and  allow  controlled  six  degree  of  freedom  hight  of  the  fuselage  is 
explored. 

II.  Single  Degree  of  Freedom  Dynamic  Model  of  RoboFly 

The  hrst  hight  of  the  Harvard  RoboFly  was  conducted  by  constraining  the  aircraft  to 
vertical  translation  on  a  pair  of  wires.  Here,  the  motion  of  the  RoboFly  fuselage  acting 
under  the  inhuence  of  a  time  varying  vertical  aerodynamic  force  that  is  a  function  of  the 
wingbeat  motion  is  modeled.  The  wings  represent  about  1%  of  the  total  vehicle  weight  and 
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Figure  1.  General  assembly  drawing  of  Harvard  RoboFly. 


Downstroke  End  of  Stroke  Upstroke 


Figure  2.  Detail  of  passive  wing  rotation  joint. 


assumed  to  be  massless  for  the  purpose  of  this  analysis.  The  equation  of  motion  for  the 
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1  DOF  Robofly  is  simply 


F^{t)  =m{x  +  g) 


(1) 


where  x  lies  along  a  unit  vector  xj  in  an  inertial  frame  that  is  taken  to  be  positive  away 
from  the  center  of  the  earth.  The  aerodynamic  force  in  the  xj  direction  is  derived  using 
blade  element  theory  for  a  triangular  shaped  wing  that  has  two  degrees  of  freedom,  namely 
angular  displacement,  0(t),  about  the  wing  root  in  the  stroke  plane,  which  is  normal  to  the 
body  x-axis  xb,  and  angular  displacement  about  the  passive  rotation  hinge  joint,  which  is 
equivalent  to  wing  angle-of-attack  a  in  still  air.  The  triangular  planform  wing  shown  in 
Figure  3  is  taken  to  be  a  rigid  flat  plate  whose  elemental  lift  at  a  spanwise  location  in  the 
local  wing  planform  plane  is  given  by 


dL  —  ^CL{a)^‘^ywpc{ywp)dywp 

(2) 

dD  =  ^CD{o)^^y‘^pc{ywp)dywp 

(3) 

where  c{ywp)  is  the  chord  at  the  spanwise  location  ywPj  which  is  a  location  on  the  wing  spar. 
The  total  lift  and  drag  on  the  wing  expressed  in  a  wing  planform  hxed  coordinate  system, 
e.g.,  RWPU,  RWPD,  LWPU,  and  LWPD  is  now  computed,  where  RWPD  is  the  right  wing 
planform  coordinate  system  for  the  downstroke,  RWPU  is  the  right  wing  planform  coordinate 
system  for  the  upstroke,  and  similarly  for  the  left  wing.  Such  a  coordinate  frame  has  an  origin 
at  the  wing  root  hinge  point  and  its  x-y  plane  is  coincident  with  the  wing  planform.  The  lift 


Figure  3.  Blade  element  computation  of  aerodynamic  forces,  moments  and  centers  of  pressure. 


is  computed  by  integrating  the  elemental  lift  over  the  span  according  to 

L=  r  dL  =  I  A  (4) 

Jo  ^ 
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Similarly,  drag  is  computed  according  to 


D=  r  dD  =  ^CD{a)^{tf  I  A  (5) 

where  I  a  is  the  area  moment  of  inertia  of  the  planform  about  the  root  and  R  is  the  length 
of  the  wing. 

Experiments  by  Sane  and  Dickenson^  for  a  dynamically  scaled  model  of  a  flapping  insect 
show  that  while  quasi-steady  estimates  fail  to  capture  the  temporal  variation  in  lift  over  a 
stroke  cycle,  it  does  capture  the  the  cycle  averaged  lift  with  reasonable  accuracy.  Best  £t 
estimates  of  quasi-steady  lift  and  drag  coefficients  derived  from  measurements  acquired  from 
180°  sweeps  of  wing  motion  as  a  function  of  angle  of  attack  are  given  by 

Cl  =  0.225  +  1.58  sin(2.13a  -  7.2) 

(6) 

Cd  =  1.92  -  1.55cos(2.04a  -  9.82) 

where  a  in  Equation  6  is  in  degrees.  For  convenience,  all  of  the  time  invariant  parameters 
are  lumped  together  according  to 

kL  =  ^CL{a)lA 

k  “  ff  i  » 

ko  =  -Cd^Ia 

Thus,  lift  and  drag  can  be  expressed  as  the  product  of  time  invariant  parameters  and  time 
varying  functions 


L  =  ki^itf 
D  =  kD(l){tf 

Note  that  the  only  variable  that  can  be  actively  manipulated  to  control  the  instantaneous 
aerodynamic  forces  is  the  angular  velocity  of  the  wing,  0(t),  and  the  forces  are  quadratic 
functions  of  this  motion  variable.  Furthermore,  because  the  vehicle  is  designed  to  mimic 
dipterian  insect  flight,  (j){t)  must  be  a  time  varying  function  that  is  equal  to  zero  at  the 
extreme  limits  of  wing  position.  It  is  assumed  that  one  can  directly  control  the  wing  position. 
This  assumption  approximates  the  physics  of  applying  a  voltage  to  an  actuator  that  induces 
a  strain  in  the  carbon  fiber  substrate  that  pushes  and  pulls  on  rigid  linkage  elements  that 
translate  the  tip  motion  of  the  bimorph  strip  into  rotational  motion  of  the  wing  root.  The 
forcing  function  that  drives  the  wing  rotation  is 

=  cos  cut  (9) 
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Assuming  that  the  frequency  of  the  oscillating  wing  is  held  constant  over  each  wingbeat 
cycle,  the  angular  velocity  of  the  wing  is  given  by 

0(t)  =  —cu  sin  cut  (10) 

Note  that  the  units  of  0(t)  are  radians  and  that  the  amplitude  of  the  wing  rotation  in  the 
stroke  plane  is  1  rad.  The  frequency  of  the  oscillator  that  drives  the  actuator  is  the  control 
input  variable.  The  proposed  control  strategy  is  based  on  the  assumption  that  the  gain 
crossover  frequency  of  the  fuselage  controller  is  much  less  than  the  trim  flapping  frequency 
required  for  hover.  If  a  non-oscillatory  control  force  were  available  and  altitude  and  altitude 
rate  measurements  were  available  for  feedback,  one  could  simply  implement  a  linear  feedback 
control  law  that  produced  the  response  of  a  damped  harmonic  oscillator.  But  since  the  system 
is  constrained  to  use  time  varying  high  frequency  oscillatory  control  inputs,  the  relationship 
between  the  time  averaged  vertical  force  and  the  control  input  is  computed.  Also  note  that 
wing  angle  of  attack  is  a  function  of  the  wing  rotation  rate  because  of  the  passive  wing 
rotation  joint  discussed  earlier. 

A.  Expression  of  Aerodynamic  Forces  in  Body  Axis  Coordinate  System 

Six  axis  systems  are  defined  to  aid  in  mapping  the  lift  and  drag  forces  acting  in  the  plane  of 
the  wing  into  body  axis  coordinates.  The  coordinate  systems  are  body,  inertial,  right  wing 
planform,  right  wing  spar,  left  wing  planform,  and  left  wing  spar.  As  shown  in  Figure  4,  a 
body  fixed  axis  system,  xb,  Ub,  zb,  is  dehned  whose  origin  is  located  at  the  center  of  gravity 
of  the  fuselage.  The  x  body  axis  lies  in  the  plane  of  symmetry  of  the  fuselage  and  the  y 
body  axis  points  out  the  right  hand  side  of  the  vehicle,  while  the  z  body  axis  points  out  the 
ventral  side  of  the  MAV.  A  right  wing  root  hxed  frame,  xrwr,  Vrwr,  zrwRi  is  dehned  that 
is  aligned  with  xbjVbjZb  but  whose  origin  is  located  at  the  right  wing  root  pivot  point.  A 
right  wing  spar  hxed  frame,  xrws,  Urws,  zrws,  is  dehned  that  rotates  through  an  angle  0(t) 
about  the  root  pivot  point.  Here,  when  0(t)  =  0,  zrws  =  —xrwr,  Vrws  is  coincident  with 
the  wing  spar,  and  zrws  completes  the  right  handed  coordinate  system.  The  transformation 
between  the  spar  and  root  axis  system  is  given  by 

Xrwr  0  0—1  xrws 

Vrwr  =  -sin0(t)  cos0(t)  0  yRws  (H) 

_zrwr\  I  cos0(t)  sin0(t)  0  J  [zrws_ 

The  right  wing  planform  frame  rotates  about  the  leading  edge  spar  by  an  angle  a  which 
corresponds  to  the  wing  angle  of  attack  when  the  vehicle  is  at  hover  in  a  quiescent  air  mass. 
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Figure  4.  Relationship  between,  inertial,  body,  root,  spar,  and  wing  planform  axis  systems. 


The  wing  plane  rotates  about  an  ideal  hinge  located  at  the  spar-root  hinge  point.  As  the 
wing  rotates  about  the  root  pivot  point,  dynamic  pressure  acting  on  the  wing  causes  it  to 
feather  into  a  direction  parallel  to  the  stroke  plane;  however,  the  magnitude  of  the  rotation 
is  limited  by  the  design  of  the  joint  as  illustrated  in  Figure  2.  When  0(t)  =  0  the  wing  hangs 
vertically,  i.e.,  a  =  90°,  under  the  influence  of  gravity  alone.  On  the  downstroke  of  the  right- 
hand  wing,  (j){t)  >  0  and  when  sufficient  dynamic  pressure  is  achieved  to  engage  the  hinge 
stop,  a  =  aiim-  Similarly,  on  the  upstroke  when  0(t)  <  0,  a  =  aum-  Note  that  a  =  aum 
on  both  the  upstroke  and  downstroke  because  positive  a  is  defined  as  a  counter  clockwise 
rotation  for  the  downstroke  and  a  clockwise  rotation  for  the  upstroke.  Figure  5  shows  the 
relationships  between  the  right  wing  spar  axis  system  and  the  upstroke  and  downstroke 
wing  fixed  coordinate  systems.  The  transformation  between  the  spar  and  downstroke  wing 
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Figure  5.  Relationship  between  right  wing  and  spar  arxis  systems  . 


coordinate  system  is  defined  as 


XrwS 

cos  a  0  sin  CK 

xrwpd 

Vrws 

= 

0  1  0 

Vrwpd 

zrws 

—  sin  a  0  cos  a 

zrwpd 

(12) 


While,  the  transformation  between  the  spar 
as 


xrws 

Vrws 

= 

_ZrwS_ 

—  cos  a 
0 


—  sin  a 


and 

0 

1 

0 


upstroke  wing  coordinate  system  is  defined 


sincK 

0 

—  cos  a 


XrwPU 

Urwpu 

_ZrwPU_ 

(13) 


With  the  relationships  between  the  body,  root,  spar,  upstroke  wing,  and  downstroke  wing 
axis  systems  established,  we  may  transform  the  instantaneous  values  of  lift  and  drag  on  the 
wing  into  the  body  axis  coordinate  frame.  Recall  that  lift  is  defined  as  the  component  of 
aerodynamic  force  perpendicular  to  the  relative  wind,  while  drag  is  defined  as  the  component 
of  aerodynamic  force  parallel  to  the  relative  wind.  If  the  air  mass  is  quiescent,  then  the 
relative  wind  is  parallel  to  the  stroke  plane  defined  by  xrws  —  Vrws-  For  fhe  upstroke,  the 
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aerodynamic  forces  expressed  in  body  axis  coordinates  are  given  by 


TpB 

Z^RWU^ 

DrwU 

Lrwu 

rpB 

Z^RWUy 

jyB  -D  RWR 

—  ^RWR^RWS 

0 

= 

—Drwu  sin(0) 

rpB 

S  RWU^. 

—LrwU_ 

_Drwu  cos(0))_ 

and  for  the  downstroke 


—Drwd 

Lrwd 

pR 

Z^RWDy 

-dR  -r  RWR 

—  ^RWR^RWS 

0 

= 

Drwd  sin(0) 

pR 

IZ^RWD^J 

_—LrwD_ 

—Drwd  cos(0)_ 

(14) 


(15) 


For  the  left  wing,  a  wing-root  hxed  frame  xlwrjUlwrj  zlwr  is  dehned  that  is  parallel  to 
xr:  Ub,  Zb  but  whose  origin  is  located  at  the  left  wing  root  pivot  point.  A  spar  hxed  frame, 
xlwSiUlwS)  zlws  is  dehned  that  rotates  through  an  angle  ^(t)  about  the  root  pivot  point. 
Here,  when  0(t)  =  0,  zrws  =  ~xlwR)  Vlws  is  coincident  with  the  wing  spar,  positive  proxi¬ 
mal  to  distal,  and  xrws  completes  the  right  handed  coordinate  system.  The  transformation 
between  the  spar  and  root  axis  system  is  given  by 


xlwr 

0  0-1 

xlws 

Ulwr 

= 

—  sin0(t)  —  cos  (f)(t)  0 

Vlws 

zlwr 

—  cos  (pit)  sin(;/)(t)  0 

Zlws 

(16) 


The  relationship  between  the  left  wing  spar  axis  and  left  wing  planform  coordinate  system 
on  the  downstroke  is 


Xlws 

—  cos  CK  0  sin  q; 

xlwpd 

Vlws 

= 

0  1  0 

Ulwpd 

_zlws_ 

—  sin  q;  0  —  cos  a 

_zlwpd_ 

(17) 


while  the  transformation  between  the  left  wing  spar  and  upstroke  left  wing  planform  coor¬ 
dinate  system  is 


Xlws 

cos  a  0  sin  a 

xlwpu 

Vlws 

= 

0  1  0 

Ulwpu 

_zlws_ 

—  sin  CK  0  cos  CK 

_zlwpu  _ 

(18) 
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For  the  upstroke,  the  aerodynamic  forces  from  the  left  wing  expressed  in  body  axis  coordi¬ 
nates  are  given  by 


rpB 

^  LWU^ 

—Dlwu 

Llwu 

pB 

^  LWUy 

-dB  -dLWR 

—  ^LWR^LWS 

0 

= 

Dlwu  sin(0) 

pB 

Flwu^. 

_—Llwu_ 

_Dlwu  cos(0)_ 

and  for  the  down  stroke 


^  LWD^ 

Dlwd 

Llwd 

^  LWDy 

■rB  T3LWR 

—  ^LWR^LWS 

0 

= 

—Dlwd  sin((;/)) 

l^LWD^J 

_—Llwd_ 

—Dlwd  cos(0)_ 

(19) 


(20) 


For  the  single  degree  of  freedom  problem,  the  cycle  averaged  force  from  both  wings, 

_ Q 

in  the  xb  direction,  ,  is  computed.  Assuming  that  the  wing  is  always  on  the  hinge 
limit  when  <j)(t)  ^  0,  then  a  =  aum  and  for  convenience  the  subscript  on  a  in  the  subsequent 
equations  is  eliminated.  By  design,  the  wings  beat  synchronously,  thus  (t>Lw(t)  =  4>Rw(t)  and 
Llwu  =  Llwd  =  Lbwd  =  Lrwu  =  L  and  these  subscripts  are  also  eliminated.  Equation  10 
indicates  that  the  wing  angular  velocity  is  cyclic  with  a  frequency  of  u  and  a  period  of  2ti /ui, 
thus  the  cycle  averaged  force  in  the  body  x-axis  direction  for  the  left  and  right  wings  is  given 
by 

(21) 


"  27r 


^LWU:,  + 


or 


^  “  m)dt 

Substituting  Equations  4  and  10  into  Equation  22  yields 


(22) 


j  r<  f  \ 
Fx  -  -l^dACiio) 


juj 


sin^  {ujt)dt 


(23) 


which  reduces  to 


=  P^ 


lACiia)  =  kLUj"" 


(24) 


Note  that  because  the  MAV  is  constrained  to  flight  in  one  direction  xj  =  xb-  Thus, 


(25) 


From  Equation  25,  the  wingbeat  frequency  that  produces  a  specihed  cycle-averaged  force  in 
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the  vertical  direction  is  given  by 


A 

- 


2F^ 

X 


pIaCl^o) 

and  a  trim  hover  frequency  can  be  computed  as 


UJq  — 


2mg 


phCda) 


(26) 


(27) 


B.  Altitude  Control  System 

If  Fxit)  could  be  manipulated  as  desired  and  if  sensors  are  available  to  measure  altitude  and 
altitude  rate,  a  linear  control  law  can  be  implemented  that  allows  the  aircraft  to  track  a 
desired  altitude  command.  A  control  law  that  produces  a  damped  second  order  response  to 
a  desired  altitude  command,  Xdes,  is  given  by 

F^{t)  =  m{-2C,ai^aX  -  ujIx  +  ulxdes  +  q)  (28) 

Of  course,  Fx{t)  cannot  be  arbitrarily  manipulated  because  it  is  constrained  to  be  periodic. 
However,  it  is  postulated  that  if  the  bandwidth  of  the  tracking  law  is  much  less  than  the 
wing  beat  frequency,  i.e.  ooa  ^  then  the  cycle  averaged  force  F^  can  be  specihed  and  used 
in  place  of  a  perfect  x-force  generator.  The  one  control  variable  that  we  have  at  our  disposal 
for  this  problem  is  the  wing  beat  frequency  uj  which  is  related  to  F^  by  Equation  26.  Figure  6 
shows  a  block  diagram  of  the  control  law.  Note  that  2^0;^  and  o;^  are  design  variables  and 
that  the  feedback  controller  manipulates  the  frequency  of  the  oscillator  that  drives  the  wing 
angular  position.  The  “Cycle  ZOH”  block  is  an  element  that  takes  the  instantaneous  value 


mg 


Figure  6.  Block  Diagram  of  Altitude  Command  Tracking  Control  System. 
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of  the  oscillator  frequency  command  signal  and  holds  it  constant  until  the  current  wingbeat 
cycle  is  complete.  Its’  purpose  is  to  preserve  the  shape  of  a  cosine  wave  over  each  wingbeat 
cycle.  Without  such  an  element,  the  control  law  can  fundamentally  change  the  shape  of  the 
wingbeat  position  wave  form  resulting  in  spikes  or  other  undesired  features. 


III.  Stability  of  Dynamic  Systems  Under  the  Influence  of 

Periodic  Control  Inputs 

In  this  section,  it  is  desired  to  analytically  determine  an  expression  for  the  altitude  of 
this  system  driven  by  a  periodic  input.  The  goal  is  to  gain  insight  into  the  problem  and  to 
determine  limitations  of  control  for  these  types  of  systems. 


A.  Analytical  Solution 

The  point  mass,  single  degree-of-freedom  equation-of-motion  for  a  flapping  wing  micro  air 
vehicle  in  hover  can  be  written  as 

x  =  -g  +  —  (29) 

m 

where  =  2L  =  2kL(j)‘^.  The  force  is  an  instantaneous  aerodynamic  force  generated 
by  the  flapping  motion  of  both  wings.  For  this  analysis,  the  angle-of-attack  of  the  wing  is 
assumed  to  be  constant  along  the  entire  span  on  both  the  upstroke  and  downstroke,  with  an 
instantaneous  change  in  angle-of-attack  at  the  end  of  each  stroke.  Thus  we  are  neglecting 
the  rotational  dynamics  at  the  end  of  each  upstroke/downstroke  where  the  wing  “flips”  in 
order  to  maintain  the  same  leading  edge,  and  thus  a  positive  angle-of-attack  with  respect  to 
the  mean  stroke  plane.  This  approximation  is  made  in  order  to  simplify  the  analysis. 

Equation  29  is  a  second-order,  ordinary  differential  equation  with  constant  coefficients. 
The  forcing  terms  are  a  step  function  with  magnitude  g  and  a  harmonic  forcing  term;  there¬ 
fore,  it  is  rather  straight-forward  to  hnd  an  analytical  solution  to  this  differential  equation. 
It  can  be  shown  that  for  the  initial  conditions  i;(0)  =  Xq  and  a;(0)  =  Xq,  the  solution  to 


x{t) 


2k  2 
—g  -\ - sm  ut 


m 


x{t)  =Xo  +  Xot  +  ^  ~ 

Inspection  of  Equation  30  suggests  that  if  xq  =  0,  then  there  is  a  unique  wing  beat  frequency 
that  will  render  the  solution  to  be  a  stable  limit  cycle  that  oscillates  about  Xq  with  an 
amplitude  kL/{Am).  The  trim  frequency  is  found  by  forcing  the  coefficient  on  the  t^  term  to 
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vanish 


(31) 


-9  kiu!^ 
2  ^  2m 


Solving  for  oo  gives  the  trim  frequency 


0 


c^o  — 


(32) 


Thus,  for  the  frequency  o^o  and  Xq  =  0,  the  vehicle  will  exhibit  a  stable  limit  cycle  described 
by 

kr 

x(t)  =  Xq  +  - —  [cos  (2ujt)  —  1]  (33) 

f  i 

with  amplitude  kiKAm)  about  the  point  Xq.  In  this  case,  for  all  frequencies  u  ^  ujq  the 
coefficient  on  will  be  non-zero;  therefore,  the  limit  cycle  will  be  unstable.  Note  that  the 
trim  frequency  calculated  in  Equation  32  is  identical  to  the  one  calculated  in  Equation  27, 
as  can  be  seen  by  recalling  that  k^  = 


IV.  Simulation  Results 

The  single  degree  of  freedom  altitude  command  tracking  system  was  simulated  in  order 
to  evaluate  system  performance.  The  vehicle  was  constrained  to  move  in  the  x  direction  only. 
The  vehicle  parameters  are  The  trim  frequency  is  equivalent  to  a  value  of  120.1742  Hz.  The 


Table  1.  Vehicle  Specifications. 


Parameter 

Value 

Units 

Mass 

60 

mg 

Width 

4 

mm 

Height 

11 

mm 

Depth 

1 

mm 

bwing 

3 

mm 

^wing 

4 

mm 

^  wing 

15 

mm 

au 

45 

deg 

OiD 

45 

deg 

COq 

755.0569 

rad 

sec. 

Ia 

935 

4 

m 

desired  altitude  command  is  a  series  of  step  functions  of  magnitude  1  m. 


Figure  7  shows 
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a  close-up  of  altitude  for  an  open-loop  simulation.  In  this  case,  the  left  and  right  wing 
frequencies  were  set  to  ujq  for  all  time.  The  vehicle  should  exhibit  a  stable  limit  cycle  given 
by  Equation  33  with  xq  =  0.  Figure  8  shows  the  altitude  from  the  simulation  as  well  as 
the  limit  cycle  from  Equation  33.  As  can  be  seen,  a  limit  cycle  does  exist  and  is  accurately 
predicted  by  Equation  33. 


Altitude  vs.  Time 


Figure  7.  Altitude  vs.  Time. 

Now,  the  feedback  loop  is  closed  and  tracking  of  a  reference  command  is  desired.  The 
reference  is  a  series  of  1  m  step  changes  in  altitude.  Figure  9  shows  the  command  and  actual 
altitudes.  Clearly,  the  system  is  capable  of  tracking  these  altitude  commands.  Recall,  the 
control  variables  are  the  fundamental  wing  beat  frequencies  of  the  left  and  right  wing.  The 
left  and  right  wing  beat  frequencies  used  to  generate  the  appropriate  altitude  are  shown  in 
Figure  10.  Relatively  large  changes  in  ujlw,  ^rw  are  necessary  in  order  to  track  the  altitude 
command.  Figures  11  and  12  show  the  altitude  tracking  and  required  wing  beat  frequency 
when  the  desired  altitude  is  a  sinusoidal  signal  of  amplitude  0.25  m  and  a  frequency  of  5 
Tracking  is  acceptable,  however,  there  is  quite  a  variation  in  wing  beat  frequency  to  generate 
a  0.25  m  amplitude  sinusoidal  wave. 
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X  10  ®  Altitude  vs.  Time  (Simulation  &  Analytical) 


Figure  8.  Altitude  vs.  Time. 

V.  Summary 

In  this  work,  an  analysis  of  the  one  degree-of-freedom  equation  of  motion  was  presented. 
An  aerodynamic  model  was  developed  using  blade  element  theory.  The  cycle-averaged  force 
in  the  inertial  x-axis  direction  (vertical)  was  computed,  from  which,  a  hover  frequency  was 
calculated.  The  point  mass  single  degree-of-freedom  equation  of  motion  was  solved  to  yield 
the  altitude  time  history.  A  stable  limit  cycle  was  achieved  with  a  suitable  choice  of  wing 
beat  frequency,  and  this  frequency  was  exactly  the  hover  frequency.  Results  showed  that  the 
solution  to  the  single  degree-of-freedom  equation  of  motion  was  identical  to  that  obtained 
from  the  full  simulation.  A  control  law,  that  produces  a  damped  second  order  response 
to  command  inputs,  was  developed.  Simulation  results  show  that  the  performance  of  the 
controller  is  acceptable. 
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Figure  9.  Altitude  vs.  Time. 
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Figure  10.  oj  vs.  Time. 
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Figure  11.  Altitude  vs.  Time. 
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Figure  12.  oj  vs.  Time. 
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